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We present an exploratory numerical study on the lattice of the color structure of the wave func¬ 
tionals of the SU(3) Yang-Mills theory in the presence of a qq static pair. In a spatial box with 
periodic boundary conditions we discuss the fact that all states contributing to the Feynman propa¬ 
gation kernel are global color singlets. We confirm this numerically by computing the correlations of 
gauge-fixed Polyakov lines with color-twisted boundary conditions in the time direction. The values 
of the lowest energies in the color singlet and octet external source sectors agree within statistical 
errors, confirming that both channels contribute to the lowest (global singlet) state of the Feynman 
kernel. We then study the case of homogeneous boundary conditions in the time direction for which 
the gauge-fixing is not needed. In this case the lowest energies extracted in the singlet external 
source sector agree with those determined with periodic boundary conditions, while in the octet 
sector the correlator is compatible with being null within our statistical errors. Therefore consis¬ 
tently only the singlet external source contribution has a non-vanishing overlap with the null-field 
wave functional. 


I. INTRODUCTION 


The knowledge of the energy levels containing infinitely heavy external sources, immersed in the quantum 
state of a pure Yang-Mills theory, is useful to understand the full QCD case |THS] • In Ref. [TO] the structure of 
the states in presence of external sources is analyzed on general grounds, within the temporal gauge, based on 
the formulation of Refs [TTT - fTTj . In particular it is shown how to define projectors on states which are irreducible 
representations of the global color symmetry, and on states where the external sources are in a given color singlet 
or octet representation. The last ones are not, in general, eigenstates of the Hamiltonian. This fact led to some 
confusion in the literature, since a state with a definite global color behavior has non-trivial projections on both 
singlet and octet external source sectors. 

It is one of the aims of this paper to clarify this issue: we show that, when the theory is quantized in 
a finite space volume with periodic boundary conditions, Gauss’s law enforces all states to be global color 
singlets. We also study in detail the properties of the ground state in presence of a static qq pair, and find that 
its wave functional projects both on the quark-singlet and on the quark-octet external source sectors. When 
homogeneous boundary conditions in the time direction are chosen, the lowest energies extracted in the singlet 
external source sector agree with those determined with periodic boundary conditions, while in the octet sector 
the correlator is consistently compatible with being null within statistical errors. 

The observables considered in this paper are correlators of Polyakov lines; it is well known that their signal 
to noise ratio is exponentially suppressed with the time-length of the lattice and with the spatial distance of 
the sources. For the pure gauge theory several techniques have been developed to overcome this problem for 
gluonic observables such as Polyakov loops, glueball interpolating operators, etc. [T5UT9] . In order to have a 
good statistical signal at large time distances, we opted for the multilevel algorithm proposed in Refs. mm- 

The paper is organized as follows. In section [TT] we collect the main aspects of the continuum theory: we 
summarize the results of Ref. m , we discuss the role of the Gauss’s law, we introduce the ” color-twisted” and 
homogeneous boundary conditions and we collect the formulas relevant for the analysis of the numerical data. 
Section m is devoted to the lattice transcriptions of the previous formulas. The results of our simulations, in 
part reported in Ref. |20j . are discussed in section IV while in section [V] we draw our conclusions. 


II. CONTINUUM THEORY 

A. External sources 

In order to compute the Feynman propagation kernel in a sector in which external q and q sources are present, 
we have to perform, in the continuum, the following operations m- 
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1. compute the propagation kernel with Aq = 0 


/*A(x,T)—A 2 (x) 

A(A 2 ,A i; T) = / Z>Aexp[-S yM (A,A 0 = 0)], (1) 

J A(x,0)=Ai (x) 

as a functional of the initial Ai and final A 2 gauge three-dimensional configurations. Eq. (jT|) describes a 
theory in which Gauss’s law is not enforced. 

2. enforce Gauss’s law in presence of qq sources, through the integration over three-dimensional gauge trans¬ 
formations n 


K(A 2 ,S 2 ,r 2 ;A 1 ,si,r 1 ]T) 


'Go 


vn n S2S1 (x,) n* ( x? -)a(a£, a i; t) , 


(2) 


where 


A^(x) = f^(x) Afj,(x) Cl(x) + i f2^(x) d^Q(x) . (3) 

In Eq. © the * denotes the complex conjugate, while si, s 2 , r\ and r 2 denote the initial and final color 
states of the q and the q sources, respectively. The domain Qq of the invariant measure integration is 
restricted to the topologically trivial gauge transformations which tend to 1 at three-dimensional infinity. 
In an infinite volume this is the condition needed in order not to eliminate from the Feynman propagation 
kernel the states with non-trivial global color content. 

It is important to notice that when the theory is defined in a finite spatial volume with periodic boundary 
conditions, Gauss’s law implies that the total color charge inside the volume carried by the gluons and the 
external sources is equal to the flux of the chromoelectric field through the boundary of the volume itself, and 
therefore vanishes identically as a consequence of periodicity. The color non-trivial states are therefore not 
present in the eigenspectrum of the finite volume theory with periodic boundary conditions, and a drastically 
different finite-size set up is needed to decide whether colored states can be excited or if their presence is 
suppressed as a consequence of the physics of color confinement. 

The kernel K in Eq. ([ 2 ]) is related to the energy eigenvalues E k and the corresponding energy eigenfunctionals 
tpk(A,s,r) by 


K(A 2 ,s 2 ,r 2 -, Ai,si,n;T) = 'Y' j e EkT ip k (A 2 , s 2 , r 2 )^(Ai, si, n) . (4) 

k 

Under a gauge transformation fl £ Qq, the states in Eq. © transform as 

i/> k (A n ,s,r) = (Cl f (x g )) aa > ( ft(xg)) r , r if> k (A,s',r ') , (5) 

and the naive scalar product of two wave functionals is infinite. It must be replaced by the Faddeev-Popov 
gauge-fixed one 


I ^MF(A)^0*(A,r,s)0(A,r,s), 

^ rs 

where 


( 6 ) 


T>M A) = A f (A)8[F(A)]8A , (7) 

and the Faddeev-Popov determinant Ap(A) is defined by (2T 

A F {A) [ VSl8[F(A a )] = 1. (8) 

Jg 0 

As a consequence of Eq. ([ 5 ]), the value of the scalar product in Eq. © is independent from the choice of the 
gauge fixing functional F(. A). In the following we assume F( A) = V • A. 
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State 

wave functional 

orbital properties 

0 at the null field conf. 

global rep 

singlet-singlet 

0 (1,1) ( A;i,j) = A) 

4>(A V ) = 0(A) 

0(0) ^ 0 

singlet 

singlet-octet 


<f> a (A v ) = <j> a (A) 

0“(O) + 0 

octet 

a-singlet 

0m*' 1 ' (A; l, j) = Sij (f>m (A) 

<j>m{A V ) = (F)0 m ((A) 

0m ( 0 ) = 0 

OL 

/3-octet 

0^’ 8) (A;i,j) = A“.^(A) 

0^A v )=< m ,(F)0^(A) 

0m ( 0 ) = 0 

7 C /3<g) 8 


TABLE I: Classification of the color states with static qq sources. Every state is labeled by the orbital and the 
source representations respectively. The symbols a and j3 indicate generic non-trivial representations. 


B. Global color transformations and classification of states 

A global color transformation V acts on the states as a unitary operator 

[U(V)^\{A ,s,r) = V ss ,V r * r ,i/j(A v ,s',r') (9) 

which commutes with the Hamiltonian |10j . If we re-organize the source indices in a matrix, the energy 
eigenstates in presence of qq sources can be written as 

^(A ,s,r) = [0(A)/ + A a 0 a (A)] sr = 0 sr (A), (10) 

where the traceless generators satisfy the normalization condition Tr[A a A b ] = S ab /2, and the color rotations in 
Eq. § can be represented as 

[W(V>](A) = ip v (A) = V^(A v )V ] = 0(A V )I + VX a V^<j> a (A v ) . (11) 

Eq. @ shows explicitly that the action of a color rotation results from a composition of a contribution coming 
from the action of V on the source indices and an “orbital” one coming from the transformation A —>• A' . 
In Tab. |T] we give the explicit classification of the color states. Several possible quantities can be computed, 
giving information about the structure of these states. In the following we classify them on the basis of the 
time boundary conditions. 


C. Time-periodic boundary conditions 


The simplest quantity to measure is the “total trace” over color source indices, i.e. the correlator of two 
Polyakov lines 


IC(R, T)= f Vvf(A) ^ K( A, s, r; A, s, r; T) = ^ d k e~ E ^ T . 

^ s,r k 


( 12 ) 


In Eq. (12 ) k runs over all the eigenstates of the Yang-Mills theory containing qq sources in all possible color 
configurations (see Table [I]); d k is the multiplicity of the fc-th level, and E k {R) is the corresponding energy 
which is a function of the relative spatial distance R = |x 9 — x^| between the two sources. The quantity 
r K(A, s, r; A, s, r; T) is gauge invariant, and in the continuum the integration over A requires a gauge 
fixing. Following Ref. m, we define the singlet and octet projectors on the color indices of the external sources 


nl 

n ; 


S2,r2,Sl,n 

8 




= 2 


A“_ A“ =d.^S r ^ - -<5, 


in terms of which we can define 


IC 1 {R,T) = 

f ©M A) 

E 

n s 2 ,r 2 , Sl ,n A(A, s 2 , r 2 ; A, si,n; T), 

J 

S2 

,T2,Si. 

,ri 

IC 8 (R,T) = 

/ ^Mf(A) 

E 

n ?2,r 2 ,»i,n^(A, s 2 , r 2 ; A, Si,ri;T). 




(13) 

(14) 


(15) 

(16) 


The integrands in Eqs. (151 and (16) are not gauge invariant and require gauge fixing. 
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D. Color-twisted boundary conditions in time 


In order to isolate the contributions from the various irreducible representations of the global color group, we 
introduce the color-twisted boundary conditions in the time direction requiring 

A(x, T) = U^A(x, 0)V , (17) 

where V belongs to the global color group SU(3). An interesting quantity is Q]]] 

K(V,R,T)= Y J VMA)V SlS2 V r * ir2 K(A v ,s 2 ,r 2 -,A,s 1 ,r 1 -,T) = Yxc*(V)e- E Z( R)T , (18) 

si,S2,ri,r2 J cx.,k 


where Xa(V) denotes the character of V in the color representation a to which the fc-th level belongs. The 
contribution from the states transforming as an irreducible representation a of the global color group are then 
given by 


K, a i y R 1 T) = dim a 


/ DVx* a (V)lC(V,R,T) 

IDV 


(19) 


where dim a is the dimension of a. Another quantity of interest is the response of the propagation kernel to an 
“orbital” color rotation, defined as 


£(v, R,T) = Y f ^Mf(A) K(A v , s, r; A, s, r; T). 

s,v ^ 


( 20 ) 


E. Homogeneous boundary conditions in time 


In Refs. 1101 13} it is observed that in presence of the homogeneous boundary conditions 

A(x,T) = A(x,0) = 0 (21) 

only the singlet-singlet and singlet-octet states survive in the kernel because they are the only ones with a non¬ 
vanishing overlap with the null-field wave functional. Under these boundary conditions correlators of Polyakov 
lines are gauge invariant, and the octet projections of the kernel on the color external state representations can 
be computed without gauge fixing. We have, in this case, 


where 


/C (0) (R, T) = /C^ 0) (R, T) + K.^ (R, T ), 


k. ( °\r,t) 

JC ( s ] (R,T) 


^(0,s,s;0,r,r;T) = Y |0fc(°)| 2 e ^ (ii)T , 


Y A(0,s,r;0,s,r;T) 


1 

3 


Y A(0,s,s;0,r,r;T) = 8^]|^(0)| 2 e ^ (fl)T , 


( 22 ) 

(23) 

(24) 


and the superscript (0) denotes homogeneous boundary conditions. 


III. LATTICE TRANSCRIPTION OF THE OBSERVABLES 


On a lattice with periodic boundary conditions we are interested in the ratio of the various kernels, as defined 
in the previous section, normalized to the analogous one without external sources. In the case of Eq. (121 the 
ratio W(R,T) is 


W(R,T) = (Tr(P(x)) Tr(P^( y ))> 5 


R - I X/y X/7 


(25) 


1 In Ref. |19| analogous boundary conditions were considered for the construction of the projectors onto irreducible representations 
of the global symmetry groups of the Yang-Mills theory. 
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where the Polyakov line P(x) is defined by 


T-l 

p(x) = n Po(x,® 0 ). 

x - 0 =0 


( 26 ) 


and (a;) is the link matrix located at the position x pointing toward the positive p direction. Under a gauge 
transformation the trace over the color index Tr(P(x)) is invariant, and, on the lattice, gauge fixing is not 
required to compute W(R,T). The correlators corresponding to the quantities in Eqs. (151 and (16) are 


lUi(P,T) = i<Tr(P(x)pt(y))) T 

Ws(R,T) = (Tr(P(x))Tr(pt(y))) r - ± <Tr(P(x)pt(y))) T , 


(27) 

(28) 


where it is understood that Coulomb’s gauge on the time slice at Xq = 0 must be fixed, see for instance Ref. [221 . 

Twisted boundary conditions in the temporal direction are introduced by requiring that 

tyx,T) =U t P M (x,0)U . (29) 

The correlators associated to the quantities K,(V,R,T) and K.(V,R,T) in Eqs. ( |T8| ) and ( |20| ) are defined as 

W(V,R,T) = (Tr(UP(x))Tr(Utpt (y)) ) T ^ , W(V,R,T) = <Tr(P(x)) Tr(pt(y))) T ^ , (30) 

where the subscript V indicates that the path integral integration is performed in presence of color-twisted 
boundary conditions, Eq. (29). 


Finally, homogeneous boundary conditions are defined by 


C/o(x,T) = [/o(x,0) = l, 


(31) 


and the definitions of the lattice correlators (R,T), w\°\R,T), (R,T) follow straightforwardly from 

Eqs. p2l), p3| and p4|. 


IV. NUMERICAL RESULTS 


The SU(3) Yang-Mills theory is discretized with the Wilson action. Numerical computations are performed 
by standard Monte Carlo techniques alternating 1 heat-bath with 5 over-relaxation steps. In the following the 
effective energy extracted from a correlator is defined as 


£(R,T) = ± In 


W(R,T- A)' 
W(R,T) \ ’ 


(32) 


where A = 4 unless T = 8 for which A = 2. 


A. Color-twisted boundary conditions 


We computed the dependence of W(V,R,T) and W(V, R,T) defined in Eq. (30), at fixed T and for different 
values of P, as a function of the octet character x 8 (V) = |Tr(U)| 2 —1. This run was carried out on a rather small 
volume, 10 3 x 4 at /3 = 6.0, in order to be able to keep under control the statistical noise without the necessity of 
using the multilevel technique. The number of generated configurations and the corresponding values of y 8 (U) 
are listed in Tab. [ITJ The values of W(V, R,T ) are independent from the character \s(U) at fixed values of R, 
see Fig. |la| This is consistent with the expectation that in presence of periodic boundary conditions every state 
contributing to the propagation kernel is a global color singlet. The states in the octet external source sector, 


which contribute to the projected Feynman kernel in Eq. (28), are therefore in an orbital-octet representation 
so to form an overall global singlet. This fact is numerically confirmed by the computation of W{V,R,T), as 
shown in Fig. [Tb] which exhibits the linear dependence on the adjoint character of V typical of states in an 
orbital octet representation. 
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Xs(V) 

# of confs. 

Xs(V) 

# of confs. 

-0.96 

10000 

5.25 

4293 

-0.46 

16000 

6.41 

2453 

0.37 

4000 

6.77 

4000 

1.12 

16000 

7.0 

4000 

2.8 

16000 

8.0 

2000 


TABLE II: Number of configurations generated and the corresponding value of the adjoint character XsiV) of 
the SU(3) matrix at the boundary. Notice that — 1 < Xs{V) < 8. 




or 

g 


Xb(V) Xb(V) 


(a) 


(b) 


FIG. 1: W(V, R, T) (a) and W(V, R, T) (b) as a function of the adjoint character y 8 (V). 


B. Periodic boundary conditions 


In order to compute the lowest energy eigenvalues of the states contributing to the correlators, the time 
extension T of the lattice has been chosen large enough for the excited states to give a negligible contribution 
with respect to the statistical errors. For this purpose, 200 independent configurations with /3 = 5.7, volume 
12 3 x T, for T = 6,8,12,16, 20 have been generated. Each of these configurations has been taken as the starting 
point for the multi-level averaging. The numerical values for the singlet and octet correlators Wi(R,T) and 
Ws(R,T) are reported in Tab |v| and Tab. VI of the Appendix Aj respectively. We computed the effective 
energies £i(R,T) and £ S (R,T) for the singlet and the octet correlator as a function of T at a given separation 
R. The results are shown in Figs. [2a| and |2b[ while their values at T = 20 are reported in Table [Till We assume 
the value of the effective energies computed at T = 20 as our best numerical estimates. 


1.6 

i | 


1.4 

i 

1 R=1 • 



R=2" 

^ 1.2 

* 

* R=3^ 

a: 





R=4a 

1.0 

R=5y 



R=6o 

0.8 

■ ■ ■ ■ 

0.6 

• • • • 


10 15 20 25 

T 


(a) 


' 

1 


1.6 

I 

f 1 


1.4 

f 

f R=1 • 


i 

R=2" 

H - 1.2 

A 

* R=3a 






R=4a 

co 

1.0 

R=5t 



R=6o 

0.8 

■ ■ ■ ■ 

0.6 

• • • • 


10 15 20 25 

T 


(b) 


FIG. 2: Effective energies £\{R,T) (a) and £&(R,T) (b) as a function of the temporal extension T of the 
lattice for R = 1 4- 6. The quantity A is equal to 4 except for T = 8 for which is 2. 
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R 

£i (11,20) 

£s(R, 20) 

1 

(5.93 ±0.02) x 10" 1 

(5.91 ±0.05) x 10" 1 

2 

(8.10 ±0.05) x lO -1 

(8.07 ±0.06) x 10" 1 

3 

1.029 ±0.009 

1.025 ±0.010 

4 

1.21 ± 0.01 

1.21 ±0.01 

5 

1.42 ±0.05 

1.42 ±0.06 


TABLE III: The values of £i(R, 20) and £s(R, 20) as a function of the distance. 


We find that the singlet-singlet and octet-octet energies coincide within statistical errors. This confirms that 
the ground state is a unique global color singlet with non-trivial components on the qq octet and singlet external 
source sectors. 


C. Homogeneous boundary conditions 


As remarked above, the use of homogeneous boundary conditions in the time direction allows the computation 
of singlet and octet source index correlators without the need of gauge-fixing. As in section |TVB[ we generated 



‘ 

I 


1.6 

t 

i i 


1.4 

T 

f R=1 • 


1 

R=2» 

K 1.2 

4 

4 R=3«. 

£ 


R-4a 

CO 



1.0 

R=5 r 



R=6o 

0.8 

■ ■ ■ ■ 

0.6 

• • • • 


10 15 20 25 

T 


(a) (b) 

FIG. 3: Effective energies £^\R, 20) (a) and £(!?, 20) (b) as a function of the temporal extension T of the 
lattice for R = 1 ± 6. The quantity A is equal to 4 except for T = 8 for which is 2. 


200 configurations on which we applied the multilevel. As compared with W(R,T), w[ 0 \R,T ) turns out to 
have a larger pre-factor as can be seen comparing the values of the kernel in Tab. |VIlT| with those in Tab. [V] In 
Figs. 


and|3b|we show the results for the effective energy £^°\r, T) of W^\R, T) and for the effective energy 
£(R,T) of the gauge invariant W(R,T) with periodic boundary conditions. The effective energies at T = 20 
are also reported in Table |TV| 

We find that the values extracted in the two cases are compatible withi n sta tistical errors, and that they 


3a 


are consistent with those computed in the previous sub-section, see Table III The correlator K^\R,T) is 
compatible with zero within statistical errors. This confirms the fact that only the singlet external source 
projection has a non-vanishing overlap with the null-field wave functional. 


V. CONCLUSIONS 

In this paper we have discussed the color structure of the Yang-Mills states in the presence of qq static quark 
sources with periodic boundary conditions in the spatial directions and periodic, color-twisted and homogeneous 
boundary conditions in the temporal one. A multilevel simulation algorithm has been necessary to obtain precise 
results on lattices corresponding to large times. From the outcome of our computations a general picture emerges 
for the structure of states in sectors with external sources. The states contributing to the Feynman propagation 















R 

£' 0) (A, 20) 

£(R, 20) 

1 

(5.90 ±0.02) x 10" 1 

(5.93 ±0.02) x 10” 1 

2 

(8.01 ±0.04) x lO -1 

(8.09 ±0.05) x 10" 1 

3 

1.011 ±0.007 

1.028 ±0.008 

4 

1.19 ±0.01 

1.21 ±0.01 

5 

1.34 ± 0.02 

1.42 ±0.04 


TABLE IV: The values of 20) and £{R 1 20) as a function of the distance. 


kernel are global color singlets only. Our interpretation is that this fact does not necessarily have a dynamical 
meaning related to confinement, but is rather due to the periodic quantization conditions imposed in the spatial 
volume. We conclude that the states are of the form referred to as singlet-singlet or octet-octet. Of course, 
any mixing between these two cases is allowed, and only dynamical considerations can determine their relative 
weight in a quantitative way. Our numerical results are consistent with the observation that only the singlet 
external source contribution has a non-vanishing overlap with the null-held wave functional. 

In view of the recent developments of the phenomenology of tetraquark (see for a review Ref. S3]) and 
pentaquark [24] states we are undertaking the study of the non perturbative effective potentials acting in these 
channels, using the techniques described in the present paper. 
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Appendix A: Numerical results for the correlators 


In this appendix we report the values of the correlators as a function of the distance between the external 
sources R and the temporal extension T computed with the multilevel algorithm. The errors are estimated 
using the jackknife method. 


R 

W 1 (R,6) 

IVl(A,8) 

IT! (A, 12) 

IT! (A, 16) 

IVl(A, 20) 

1 

(2.666 ± 0.007) x 10~ 2 

(8.19 ±0.03) x 10" 3 

(7.72 ±0.03) x 10" 4 

(7.25 ± 0.03) x 10" 5 

(6.82 ±0.03) x 10“ 6 

2 

(6.28 ±0.04) x 10 -3 

(1.254 ±0.009) x 10~ 3 

(5.05 ± 0.04) x 10" 5 

(2.02 ±0.02) x 10 -6 

(8.04 ±0.09) x 10“ 8 

3 

(1.42 ±0.02) x 10 -3 

(1.82 ±0.02) x 10 -4 

(3.11 ± 0.04) x 10 -6 

(5.28 ±0.09) x 10" 8 

(8.86 ±0.18) x 10" 10 

4 

(3.78 ±0.06) x 10" 4 

(3.29 ±0.06) x 10" 5 

(2.68 ±0.06) x 10" 7 

(2.19 ± 0.06) x 10 -9 

(1.79 ±0.06) x 10 -11 

5 

(1.15 ±0.04) x 10 -4 

(6.6 ±0.2) x 10" 6 

(2.7 ± 0.1) x 10“ 8 

(1.03 ±0.07) x 10” 10 

(4.2 ±0.5) x 10 -13 

6 

(6.3 ±0.4) x 10" 5 

(2.5 ±0.2) x 10 -6 

(6.0 ±0.9) x 10“ 9 

(3.7 ±6.9) x 10" 12 

(5.1 ±4.9) x 10~ 14 


TABLE V: Singlet projected correlator computed with the gauge-fixing and periodic boundary conditions. 


R 

W a {R, 6) 

IMA, 8) 

IMA, 12) 

IMA, 16) 

IMA, 20) 

1 

(2.42 ± 

0.02) x 10" 3 

(7.11 ± 

0.06) 

x 10” 4 

(6.63 

± 

0.06) x 10" 5 

(6.22 ±0.06) 

x 10 -6 

(5.88 ±0.06) x 10“ 7 

2 

(1.92 ± 

0.02) x 10" 3 

(3.64 ± 

0.03) 

x 10~ 4 

(1.44 

± 

0.01) x 10" 5 

(5.74 ±0.07) 

x 10 -7 

(2.31 ±0.03) x 10“ 8 

3 

(9.4 ± 

0.1) x 10“ 4 

(1.09 ± 

0.01) 

x 10 -4 

(1.82 

± 

0.02) x 10~ 6 

(3.04 ±0.05) 

x 10 -8 

(5.18 ±0.10) x 10 _1 ° 

4 

(4.8 ± 

0.1) x 10“ 4 

(3.56 ± 

0.07) 

x 10 -5 

(2.76 

± 

0.05) x 10" 7 

(2.21 ±0.06) 

x 10 -9 

(1.81 ±0.05) x 10' 11 

5 

(2.7 ± 

0.2) x 10“ 4 

(1.27 ± 

0.04) 

x 10~ 6 

(4.6 

± 

0.2) x 10 -8 

(1.9 ±0.2) x 

io- 10 

(7.2 ±0.1) x 10 -13 

6 

(2.1 ± 

0.1) x 10“ 4 

(7.3 ± 

0.4) x 10' 6 

(1.4 

± 

0.3) x 10" 8 

(2.2 ± 1.9) x 

10" 11 

- 


TABLE VI: Octet projected correlator computed with the gauge-fixing and periodic boundary conditions. 
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R 

W(R, 6) 

IT (A, 8) 

W(R, 12) 

W{R, 16) 

W(R, 20) 

1 

(2.908 ± 

0.006) x 10" 2 

(8.90 

± 

CO 

1 

O 

X 

c? 

o 

o 

(8.39 

± 

0.03) x 10~ 4 

(7.87 ±0.03) x 10“ 5 

(7.41 ± 0.04) x 10 -6 

2 

(8.20 ± 

0.04) x 10 -3 

(1.62 

± 

0.01) x 10 -3 

(6.50 

± 

0.04) x io- 5 

(2.59 ± 0.02) x 10“ 6 

(1.04 ± 0.01) x 10 -7 

3 

(2.36 ± 

0.02) x 10 -3 

(2.92 

± 

0.03) x 10“ 4 

(4.94 

± 

0.06) x 10^ 6 

(8.32 ±0.13) x 10“ 8 

(1.40 ± 0.03) x 10 -9 

4 

(8.6 ± 

0.2) x 10 -4 

(6.90 

± 

0.11) x 10 -5 

(5.43 

± 

0.09) x 10“ 7 

(4.4 ± 0.1) x 10“ 9 

(3.59 ±0.10) x 10~ u 

5 

(3.9 ± 

0.1) x 10“ 4 

(1.93 

± 

0.06) x 10“ 5 

(7.3 

± 

0.3) x 10' 8 

(3.0 ±0.2) x 10 _1 ° 

(1.1 ±0.1) x 10 -12 

6 

(2.8 ± 

0.1) x 10 -4 

(9.8 

± 

0.5) x 10~ 6 

(2.1 

± 

0.3) x 10' 8 

(2.6 ±2.0) x 10 -11 

- 


TABLE VII: Total trace correlator in presence of periodic boundary conditions. 


R 

IV^ (A, 6) 

W[ 0 \R, 8) 

W^ 0 \R, 12) 

W[ 0) (R , 16) 

W 1 (0) (R, 20) 

1 

(5.065 ± 0.003) x 10~ 2 

(1.154 ±0.003) x 10~ 2 

(1.424 ±0.004) x 10~ 3 

(1.346 ± 0.005) x 10“ 4 

(1.265 ± 0.005) x 10~ 5 

2 

(2.713 ±0.003) x 10~ 2 

(5.17 ±0.03) x 10" 3 

(1.96 ±0.01) x 10" 4 

(7.86 ±0.06) x 10" 6 

(3.15 ±0.03) x 10" 7 

3 

(1.684 ± 0.003) x 10~ 2 

(1.97 ±0.02) x 10“ 3 

(2.91 ±0.03) x 10" 5 

(4.89 ±0.07) x 10" 7 

(8.38 ±0.13) x 10" 9 

4 

(1.322 ±0.003) x 10~ 2 

(1.04 ±0.02) x 10 -3 

(6.33 ±0.11) x 10" 6 

(5.01 ±0.10) x 10" 8 

(4.19 ±0.10) x 10 _1 ° 

5 

(1.172 ±0.003) x 10~ 2 

(6.87 ±0.16) x 10" 4 

(1.70 ±0.04) x 10" 6 

(6.4 ± 0.2) x 10 -9 

(2.8 ±0.1) x 10 -11 

6 

(1.130 ±0.003) x 10~ 2 

(5.93 ±0.16) x 10" 4 

(8.5 ±0.3) x 10“ 7 

(1.7 ± 0.1) x 10“ 9 

(5.1 ±0.8) x 10- 12 


TABLE VIII: Singlet projected correlator computed in presence of homogeneous boundary conditions. 
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